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limit the intensity of the fermion field coincides with the transition probability in the two neutrinos 
system interacting with moving and polarized matter. 
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I. INTRODUCTION 

The description of the mixed fermions evolution at- 
tracts considerable attention after the experimental con- 
firmation of solar neutrino oscillations 0,12 • The major- 
ity of the neutrino oscillations studies involve the quan- 
tum mechanical approach to the description of the neu- 
trino wave function evolution (see especially Ref. 0). 
Despite of the fact that quantum mechanics allows one to 
establish the main properties of the neutrino oscillations 
process this method of the neutrino oscillations treatment 
has several disadvantages. Neutrinos are usually sup- 
posed to be scalar particles without reference to the mul- 
ticomponent single neutrino wave function. The famous 
Pontecorvo formula (see Ref. Q), which is in use in many 
theoretical and experimental studies of neutrino oscilla- 
tions, is valid only for ultrarelativistic particles. How- 
ever, the exact theory of the considered process must be 
applicable for neutrinos with arbitrary energies. Quan- 
tum mechanical approach also does not make it clear if 
mass or flavor eigenstates bear more physical meaning. 
Therefore one can see that a theoretical model of neutrino 
oscillations, which would overcome the mentioned above 
difficulties, should be put forward. There were numer- 
ous attempts to construct the appropriate formalism for 
neutrino flavor oscillations in vacuum. The quantum field 
theory was applied to this problem in Refs. d, IE @ • 
The authors of these papers reproduced the Pontecorvo 
formula and discussed the corrections to this expression. 
Recently we revealed in Ref. that neutrino flavor os- 
cillations in vacuum could be explained in frames of the 
classical field theory. 

It was also realized that neutrino interactions with ex- 
ternal fields can drastically change the picture of the 
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oscillations process. For example, it was discovered in 
Refs. [jjj, uM that the transition probability can achieve 
great values if a neutrino interacts with background mat- 
ter by means of weak currents. Thus we should develop 
now not only the appropriate theory of neutrino oscil- 
lations in vacuum, but also include in out treatment 
possible effects of neutrino interactions with external 
fields. During the last three decades the approaches for 
the theoretical substantiation of the Mikheyev-Smirnov- 
Wolfenstein (MSW) effect were developed. Among them 
we can distinguish Refs. 0,^1 i n which the neutrino rel- 
ativistic wave equations in dense matter were analyzed. 
The S'-matrix approach was used Ref. 14] to account for 
the MSW effect. The influence of moving and polarized 
matter was described in Refs. [l5llT^ |. 

The purpose of the present work is to provide deeper 
understanding of the neutrino flavor oscillations phe- 
nomenon. The approach developed in our paper can not 
only reproduce the Pontecorvo formula for the transition 
probability, but also give clear physical explanation of 
the corrections to this expression, which are widely dis- 
cussed now (see, e.g., review [T^ | and references therein). 
The analysis used in this article is based on classical field 
theory methods. We study the evolution of coupled clas- 
sical fermions under the influence of external axial- vector 
fields. A classical fermion is regarded as a first quantized 
field because the Dirac equation, 

1 dtp _ , imc - , 

+ aVt/j+ — Pip = 0, 
c at n 

in which we use the common notations for the gamma 
matrices a = 7°7 and (3 = 7 , does contain Plank con- 
stant H, however the wave function tp is supposed to be 
the non-operator object in our approach. Therefore we 
do not involve the second quantization in the present 
work. Note that the discussion of the first quantized 
neutrino fields is also presented in Ref. [Tsj . 

In Sec.[n]we start from the flavor neutrino Lagrangian 
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which accounts for the interaction with external axial- 
vector fields. Then we derive the basic integro-differential 
equations for the "mass eigenstates" which exactly take 
into account both Lorentz invariance and the interaction 
with external fields. These equations are also valid in 
(3+1) dimensions. The perturbation theory is used for 
the analysis of the obtained equations. In Sec. IHII we 
discuss the neutrino fields distributions and obtain the 
zero order term in their expansion over the external fields 
strength. This result corresponds to vacuum neutrino os- 
cillations. In Sec. II VI we get the first order correction to 
the neutrino field intensity in vacuum and show that our 
formula is identical to the expression for the neutrino 
transition probability obtained in the quantum mechan- 
ical approach for ultrarelativistic neutrinos. This case 
corresponds to neutrino flavor oscillations in moving and 
polarized matter. Finally we discuss our results in Sec.lVl 



II. GENERAL FORMALISM 

Without losing generality we can discuss the evolu- 
tion of the two coupled fermions system (1^1,^2)- These 
fermions are taken to interact with the external axial- 
vector fields fi 2 - The Lagrangian for this system is ex- 
pressed in the following form 



fc=l,2 
k=l,2 



where Gf is the Fermi constant and 



Ik ' 



(2.1) 



where g is the coupling constant (g* is the complex con- 
jugate value), 7^ = 7^(1 + 7 5 )/2 ( 7 5 = -WlV), 
and 

£o(vk) = ^fe(«7 M ^ - ™-k)Vk- 

is the free fermion Lagrangian, rttfc are the masses of the 
corresponding fermions Vk- 

One of the possible examples of the fermions up. is the 
system of neutrinos belonging to different flavor states. 
In this case we can identify the first fermion in Eq. (|2.1|) 
with a muon neutrino or a r-neutrino v r and the sec- 
ond one with an electron neutrino v e . These neutrino 
types are known to interact with matter composed of 
electrons, protons and neutrons by means of the elec- 
troweak interactions. Note that an electron neutrino in- 
teracts with background fermions via both charged and 
neutral weak currents whereas a muon or a r-neutrino 
are involved only in the interaction through weak neu- 
tral currents. Thus the external axial-vector fields f^ 2 
can be expressed in terms of the hydrodynamical cur- 
rents jy? and the polarizations of different fermions in 
matter (see, e.g., Refs. \1&L fl9t 120] ) . 
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Here is the third isospin component of the matter 
fermion /, is its electric charge and 9w is the Wein- 
berg angle. The hydrodynamical currents and the polar- 
izations are related to the fermions velocities Vf and the 
spin vectors £f by means of the following formulas 



3f = (nf,n f v f ), 



(2.4) 



X f = ( «/(C/V/),n/C/ A /l 
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The detailed derivation of Eqs. I|2. 211 - 12. 4|l is presented in 
Refs. HHllll^. 

Following the results of our previous work we will 
study the evolution of the system (|2.1|) by solving the 
Cauchy problem. Let us choose the initial conditions in 
the form 



i/i(r,0) = 0, i*(r,0)=£(r), 



(2.5) 



where £(r) is the known function. If one considers the 
fermions as flavor neutrinos, the initial conditions in 
Eq. (|2.5|) correspond to the common situation in a neu- 
trino oscillations experiment, i.e. v^ T are absent initially 
and v e has some known field distribution. We will be in- 
terested in searching for the fields distributions Vj- (r, t) 
for t > 0. 

In order to solve the Cauchy problem we introduce the 
new set of the fermions (V'i, ^2) by means of the matrix 
transformation, 



cos 9 sin 9 \ ( ipi 
-sin0 cos 9 J 1-02 



(2.6) 



The mixing matrix (which is parameterized with help of 
the one angle 9) in Eq. I|2.6|l is chosen in the way to 
eliminate the second and the third terms in Eq. (|2.1|l . 
If we had studied the evolution of our system without 
the external fields fi 2 , i-e. in vacuum, the fermions ipk 
would have been called mass eigenstates because they 
would have diagonalized the Lagrangian and therefore ipk 
would have had definite masses. In our case (/(* 2 7^ 0), if 
we simplify the vacuum mixing terms, the matter term 
becomes more complicated compared to Eq. (|2.1f) . The 
Lagrangian rewritten in terms of the fields ipk can be 



expressed in the following way, 



fe=l,2 



where 



+ ^27>2(c 2 / 2 AI + S 2 /f) 

+ scfa-yflto + ^27^i)(/r - #)], (2.7) 



It should be noted that the masses m k of the fermions 
ip k are related to the masses xtv k by the formula, 



mi 



c 2 mi 



s 2 m 2 , 



m 2 



c 2 m 2 



s 2 mi, 



In Eq. 1)2. 7|l we use the notations c = cos 9 and s = sin#. 

Eq. 1)2.7)1 has some advantages in comparison with 
Eq. 1)2.1)1 despite of the more complicated matter inter- 
action term. The terms gD 2 v\ and g*v\V2 in Eq. 1)2. 1)1 . 
which are responsible for the vacuum oscillations cannot 
be treated within the perturbation theory. In order to 
describe the flavor changing processes one has to take 
into account these terms exactly. In Eq. I|2.7|l we have 
the additional interaction terms which can be analyzed 
in the common way with help of the perturbation the- 
ory if the strength of the fields /f 2 is supposed to be 
weak. The criterion of the external fields weakness is 
discussed in details in Sec. II VI It is also possible to as- 
sume the convenient time dependence (sometimes it is 
necessary to consider the adiabatic switching-on of the 
interaction) of the fields /f 2 in Eq. (|2.7|l . On the con- 
trary we cannot " switch-on" or " switch-off" the constant 
g in Eq. 1)2.1)1 at some moments of time because it is re- 
lated to the properties of a theoretical scheme underlying 
the phcnomcnological model of neutrino oscillations. 

Now let us discuss the evolution of the system ()2.7|) 
with the initial conditions 

^i(r,0) = -s£(r), ^(r, 0) = c£(r), 

which follow from Eqs. 1)2.5(1 and 1)2.6(1 . We also suppose 
that the external fields /(* 2 are weak and one is able to 
take them into account in the lowest order of the per- 
turbation theory. If we had solved this problem using 
the quantum field theory, we would have calculated the 
contributions of the four Feynman diagrams shown on 

Fig.m 

We can always rewrite the Dirac equations which result 
from Eq. 12.7(1 in the form, 



iip 1 = (H 1 +V 1 )ip 1 + Vip 2 , 
iip 2 = {H 2 + V 2 )ip2 + V<if> 1 , 



(2. 



where V 1 = /3 7 £(c 2 /f + s 2 tf), V 2 = Prf(c 2 tf + S 2 /f), 
V = sc/?7^(/C — $2) an d Hk = —iaS7 + l3m k are the free 
fields Hamiltonians. We are searching for the solutions 
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FIG. 1: Feynman diagrams contributing to the interaction of 
the fermions ip k with the external axial- vector fields /f 2 . 



of Eqs. 1)2.8)1 in the following way, 



V>fc(r,i) 



d 3 p 

(2tt) 3 / 2 ' 



MP>*)*£ P (*)+&k(P>W. D (s)]> (2-9) 



where *£ p (a:) = u k (p)e~ ipkX and % tP (x) = v k {p)e lpkX 
are the basis spinors, = (t,r), — (£ k ,p) and 

£k = \Jp 2 + m\ is the energy of the fermion ipk. The 
coefficients etfc(p,i) and b k (p,t) are not the creation and 
annihilation operators since we are using here the classi- 
cal field theory. The values of these functions should be 
chosen in the way to satisfy the initial conditions 1)2.5)1 . 
Note that there is the additional time dependence of these 
functions in contrast to the case of the flavor changing 
process in vacuum (see Ref. @). 

Using the orthonormality condition of the basis spinors 
ty+ p (x) = u k {p)e- ipkX and ^ p {x) = v k (p)e lPkX we ob- 
tain from Eqs. 1)2.8)1 the new integro-differential equations 
for the functions a k (p,t) and b k (p,t), 

»*i(P>*) =(2^372 (/ d 3 rHf+1,(x)ViM*,t) 
+ J d 3 r^ p (x)VMr,t?j, 



ia 2 (p,t) 



(2tt) 3 / 2 

+ J d 3 v*-J p (x)VMr,t)^, 

d 3 r^ p {x)V 2 Mr,t) 



(2tt) 3 / 2 
+ J d 3 r*++(a)V^i(r,«)Y 



4 



^(p, *) =J2^ { J d * T *lU x ) V *M*, t) 

+ j d 3 r^ p (x)VM^t)\ (2-10) 

These equations correctly take into account both the 
Lorentz invariance and the interaction with the exter- 
nal fields /f 2 . However if we suppose that these external 
fields are rather weak, it is possible to look for the solu- 
tions of Eqs. (|2.1U|) in the form of the series, 



a k (p,t) =af(p)+a£\p,t) + .. 



b k (p,t) =&i 0) (p) + 6n P ,i) + 



(2.11) 



Eqs. H2.11f) mean that the fields distributions can be also 
presented in the form of the series ^(r, t) — ip^\r, t) + 
tp^ (r, i) + • ■ • . Note that the coefficients af* 1 (p) and 



^(p); which correspond to the function ip£'(r,t), do 

not depend on time. The functions ipf^ ( r ' are respon- 
sible for the evolution of the considered system in vac- 
uum, i.e. at /f 2 = 0. These functions in two dimen- 
sional space-time have been found in our previous work 
where the analogous Cauchy problem has been solved 
in the explicit form in (l+l)-dimensional case. However 
to describe the evolution of our system with the non-zero 
external fields in the (3+l)-dimensional space-time we 
should study the vacuum case also in (3+1) dimensions. 



(0), 



where 

S k (-p,t) 
and 



cos£ k t - i Sm f ht (ap + (3m k ) (0), (3.3) 

Ok 



Up) = J d 3 pe- ipr £(r), 

are the Fourier transforms of the functions Sk{r,t) and 

Now let us choose the initial condition. We suppose 
that the initial field distribution of v 2 is the plane wave, 
i.e. £(r) = e lujr ^Q, where £o is the normalization spinor. 
The Fourier transform of this function can be simply 
computed, f(p) = (2tt) 3 <5 3 (u; - p)£ . Using Eqs. fFfy 
we get the fields distributions in the (3+l)-dimensional 
space-time for the plane wave initial condition, 



cos[£i(w)£] 



. sinffi (u))t] , . . 



cos[£ 2 (uj)t} 



i — ——^ — (au + (3m 2 ) 



£o, (3.4) 



III. EVOLUTION OF THE SYSTEM IN 
VACUUM 

To study the behavior of ^4 (r, t), i.e. the evolution 
of our system in vacuum, we use the results of Ref. |9j 
where it was revealed that the fields distributions for the 
given initial conditions have the form, 

i>{°\r,t)=-s J d 3 r%(r'-r,t)H/3)£(r'), 

4°\r,t) =c J d 3 r%(r' - r, t)(-i/?)£(r'), (3.1) 

where <S/c(r,t) is the Pauli- Jordan function for the 
fermion ipf. (see, e.g., Ref. |5l|). It should be noted that 
Eqs. H3.1fl are the most general ones and valid in (3+1) 
dimensions. Contrary to the approach of Ref. here 
we use the momentum representation because the inte- 
grations in Eqs. I|3.1|) are rather cumbersome in the co- 
ordinate representation. Thus one rewrites these expres- 
sions using the Fourier transform of the initial conditions. 
Eqs. H3.1(l now take the form, 

4°\r,t)=-s J ^e^-p^H/^p), 
4°\r,t) = c J e ^5 2 (-p,t)(-i/3)e(p), (3.2) 



where £k{^>) — v/w 2 + writ. 

In the following we discuss the case of the rapidly os- 
cillating initial conditions, i.e. lj 3> rn\ t 2. One obtains 
from Eqs. I|3.4|l the fields distributions for to ^> mi,2 in 
the following form, 



4 0) M) 



se %MT {eas\£ x {u)t\ 

- i(an) sin[£i(a;)i])£o, 
4°\r,t) =ce lujr (cos[£ 2 (u;)t} 

- i(an) sm[£ 2 (ui)t])£o, 



(3.5) 



where n = u>/u> is the unit vector in the direction of the 
initial field distribution momentum. The fermion v\ is 
absent at t = 0. Therefore it would be interesting to 
examine the field distribution vf^ (r, t) at t > 0. Using 
Eqs. (|2.6|) and l|3.5|) we obtain, 

^ 0) (r,i) =c4 0) + s4 0) = sin26>sin[A(w)i] (3.6) 
x {sin[cr(w)<] +i(an)cos[cr(w)t]}e lwr $ , 



where 



a(oj) — > w H 



4lJ 



AH 



£i(w) — ^2(w) m 2 



ml 



4w 



Am 2 
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The measurable quantity of a classical spinor field is the 
intensity. With help of Eq. I|3.5I) one gets the intensity of 

the fermion in the following form, 

l[°\t) =\v[ 0) (r,t)\ 2 =sin 2 (20) sin 2 [AHi] 

x ^| sin[cr(w)t] + i(an) cos[er(w)t]| 2 <!;o 

= sm 2 (2#)sin 2 ( ^-t) , (3.7) 



which reproduces the Pontecorvo formula in the (3+1)- 
dimensional space-time since we can regard the intensity 
of the fermion V\ as the transition probability in two 
neutrino system. Eq. (|3.7I) also generalizes the result of 
our previous work 9] where the analogous expression was 
derived in (1+1) dimensions. 



I 

IV. INTERACTION OF THE SYSTEM WITH AN EXTERNAL FIELD 



In order to proceed in our studying of the two neutrino system evolution under the influence of the external fields 
f^ 2 we discuss the further correction to the vacuum case. The first order corrections to Eqs. i|3.2|) can be derived 
from Eqs. (|2.9() and (|2.10l) and have the form, 



d 3 p 



= ipr 



(2tt) 3 2£ 



-{£i[s^i(6+e- !;£l * + 6r e +l£l *) - cV{&+ 2 e~ t£lt + &^ 2 e +l£lt )] 



+ (ap + Pm 1 )[sV 1 (e+ 



,-iEit 



+i£lt\ 



4 1 \r,t)=-i 



d 3 p e 



ipr 



-{£ 2 [cV 2 (6+e- l£2t + & 2 e +l£2t ) - sV{e+ 1 e- l£2t + G 21 e +l£2t )] 



where 



(2tt)3 2£ 2 

+ (ap + (3m 2 )[cV 2 (&+e- l£2t - & 2 e +l£2t ) ~ sV(6+e~ 4£2t - 6 21 e + ' l£2t )}}(-i(3)£(p), 

©12=/ e ±l£lt S 2 (-p,t)dt, &fi= f e ±i£ ^S x {-p,t)dt, ef <2 = [ e ±l£l2 *5 li2 (-p,t)dt. 
Jo Jo Jo 



(4.1) 
(4.2) 



When we derive Eqs. (|4.1(l we suppose that a^(p,0) = 6i^(p,Q) = 0. It means that at t = the fields distributions 
are determined by Eqs. (|3.2() . One can find the explicit form of the functions given in Eqs. I|4. 2J1 using Eq. I|3.3(l . 



1 L±«rt sin ^ , „±iAt siaAt ap + (3m 2 ( ±e±i0 . t sin at ^ ^ ±iAt sin At 
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®f,2 = 







+ e 



A £ 2 
TiAt sinAt ap + f3rrii 



Si 



±e 
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sin at 



+ e J 



a 



+ e 



=piAt 



A 

sin At 



±2i£ li2 t 



+ t]+i 



ap + (3m h2 f 1 ±2i£l 

Sl,2 



where 



Sl+S2 A 

O = — - — , A 



Si-S 2 



2 ' + t) }(*/?), 



(4.3) 



On the basis of Eqs. i|4.1|) and (|4.3|l we obtain the expressions for the first order corrections to the vacuum case in 
the following form, 



^\v,t) =i 



d 3 p e lpr 
(2tt) 3 2£^ 
sin at 



£1 



sVi 



sin£it .ap + /3mi . 

rcoscit — j tsiniit 



Si 



£1 



cV 



sin At ap + [3m 2 
cos At H — cos at + i sin at sin At 



+ (ap + /3mi) 
-cV 



sVi < —it sin fit 



£2 

ap + /3mi 
£1 



1 1 

a ~ A 



tcos £it 



sin fit 
£1 



4 1) (r,t)=-z 



— i sin at sin At 
d 3 p e 4p 



1 1\ ap + (3m 2 /sin At sin at 

— — cos at cos At 



(2tt) 3 2£ 2 



£2 



cV: 



a A 

$in£ 2 t 
£2 



ap + 0m 2 

tcos£ 2 t — 1 1 sm£ 2 t 

£2 
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sV 



sincrt . sin At ap + /?mi 
cos At H : — cos at — i sin at sin At 



A 



£1 



l l 

a + A 



+ (ap + (3m 2 ) 



cV 2 



-irsint2iH ^ t cos e 2 t ^ — 



— sU s « sin at sin Af ( — 

a A 



1 l\ ap + Qm\ /sinAi sin at 

— — cos at cos At 



(4.4) 



Note that these expressions are valid for the arbitrary initial conditions £(p) and exactly take into account the 
Lorentz invariance. It should be also mentioned that along with the harmonic functions there are several terms in 
the integrands which linearly depend on time. Therefore for Eqs. H4.4|l to be meaning-bearing one has to assume that 
the potentials Vi t2 and V are rather weak, i.e. we study the interaction of our system with weak external fields (see 
Eq. flOj l below). 

The integrations in Eqs. (|4.4|l are rather complicated for the arbitrary initial conditions. That is why we choose 
the function £(p) analogously to Sec. IHII i.e. we again suppose that £(p) = (2it) 3 5 3 (uj — p)£o- The integrations over 
momenta are eliminated with help of the <5-functions. We also consider high frequency approximation, i.e. u) 3> mi^. 
Finally we obtain the following expressions for the fields distributions of the fermions ipk, 

4 1} (r,t) =e*"i 



st[Vi + (an)Vi(an)](cos[£iH*] — i(an) sin[£i(w)t]) 
sin[A(w)t] 



AH 



[V + (an)V(an)](cos[a(uj)t} — i(otn) sin[a(uj)t]) >£ , 



V^ 2) (r,<) = - e^ r -|ct[y 2 + (an)V 2 (an)](cos[£ 2 (uj)t\ - t(an) sin[£ 2 H*D 

- s Sm ^^ [V + (an)V(an)](coB[<7(tt)t] - i(an) sin[a(w)<]))eo- 



(4.5) 



On the basis of Eqs. i|2.6[l and l|4.5|l we can derive the first order correction to the field distribution of the fermion v\ 
in the form, 



vi\r,t) =-sin26>e"" r -<J cos 2(9 



sin A(ui)t 



(Fi — F 2 )(cos[a(u))t} — i(an) sm[a(u!)t]) 



t[Fi{c 2 (cos[fi(w)t] - i(an) sin[£i(w)t]) - s 2 (cos[£ 2 (ui)t] - i(an) sin[£ 2 H*])} 
F 2 {c 2 (cos[£ 2 (uj)t} - i(an) sm[£ 2 (u)t\) - s 2 (cos[£iHi] - i(an) sm[£x(u))t])}] W (4.6) 

I 



where Fi, 2 = [fl 2 - (fi, 2 n)(En)](l + 7 5 ) and £ = - 7 5 a. 

To calculate the intensity of the field v\ one should take 
into account that the final expression for the intensity 
must contain only terms linear in external fields. There- 
fore the first order correction to the intensity should be 
calculated with help of the formula 



,(o)t„(i) 



,(i)t„(o) 



Using Eqs. and l|4.6|) we get the expression for i" l7 

lP(t) = sin 2 (2(9)cos26>sin[A(w)t] (4.7) 
1 /sin[A(cj)i] 



AH 



tcos[A(w)f] 



x (([/"(an) - (f 2 n)] - [/°(an) 



(fm)])(l+7 5 )>- 



In Eq. (|4.7|) we use the notation ((...)) = £o(- ••)&)• 
To compute the mean value with help of the nor- 
malization spinor £o we can suppose that £(r) = 



exp(— iE„ 2 t)^(r)\ t ^,o- Then we notice that for spinors 
corresponding to high energies one has the obvious iden- 
tities (1 +7 5 )£o ~ 2£ and (J(an)( « 1. Putting to- 
gether Eqs. (|3.7(l and (|4.7|l we obtain the final expression 
for the intensity of the fermion i>\, 

h(t) =jW(t)+jW(t) = 

sin 2 (26>) j sin 2 [A(w)£] + cos 29 sin[A(w)t] 

' sin[A(uj)t\ 



t cos [AH*] 



x([/ 2 °-(f 2 n)]-[/ 1 °-(f 1 n)]) 



(4.8) 



Using Eq. (|4.8|l it is possible to define the scope of the 
applied method, i.e. we can evaluate the strength of ex- 
ternal fields necessary for the perturbative approach to 
be valid. With help of Eq. (|4.8|l one obtains the inequal- 
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ities, 



Acos20 < A(w), Aicos26 ) <l, 



(4.9) 



where A = [/ 2 ° - (f 2 n)] - [/£ - (fin)]. If Eq. g2J is 
satisfied, the contribution of external axial-vector fields 
to neutrino flavor oscillations is small compared to the 
vacuum term. It should be noted that Eq. (|4.9(l is valid in 
the ultrarelativistic case. For neutrinos with £&(w) ~ 
we should rely on Eqs. i|4.4[l rather than on Eqs. I|4.5|) . In 
this case the condition of our method applicability will 
be different from Eq. (|4.9|l . 

Now let us compare Eq. I|4.8|l with the neutrino transi- 
tion probability formula. Flavor neutrinos are considered 
to interact with external axial-vector fields as it is de- 
scribed in Eqs. (|2.2[> - (|2.4[1 . Then the probability to find 
muon or T-neutrinos in the electron neutrinos beam in 
presence of moving and polarized matter is expressed in 
the following way (see, e.g., Refs. 003) 



irt 

P» e ^v„ r (t) = sm 2 (26W) sin" ( — 

Lett 



(4.10) 



where 



sin 2 (26> off ) 



A 2 (w) sin 2 (2(9) 



[A(w) cos 26» - A/2} 2 + A 2 H sin 2 (26») ' 
is the definition of the effective mixing angle and 



(4.11) 



— = J [A(w) cos26> - A/2} 2 + A 2 H sin 2 (26), (4.12) 
Lett v 

is the definition of the effective oscillation length. 

To discuss the weak external field limit in Eqs. (|4.10|) - 
(|4.12l) we should expand Eqs. (|4.11|) and l|4.12|l over small 
parameter A [see Eq. I|4.9|) ]. As a result one gets 



sin"(26> cff ) « sin" (261) 1 



Acos26> 
tt/L cS wA(cj) - (A cos 260/2. 



(4.13) 



It is also necessary to expand the time dependent factor 
in Ea. |@Znfl) , 



sin 



;sin 2 [2A(w)t] 

- At cos 29 sm[A(uj)t] cos[A(w)t]. (4.14) 

Note that Eq. (|4.14|l is valid while At cos 29 < 1, that co- 
incides with the second inequality in Eq. I|4.9|l . With help 
of Eqs. 1)4. 13[) and (|4.14l) the neutrino transition proba- 
bility is reduced to the form, 

Pu^u^.M =sin 2 (26»)^ sin 2 [AH£] + cos 29 sin[A(w)i] 

'sin[A(cj)tl ... . . 

1 v ; J -tcos[A(u)t] 



x([/ 2 °-(f 2 n)]-[/ 1 °-(f 1 n)]) 



which coincides with Eq. (|4.8() . This comparison shows 
that neutrino flavor oscillations in weak axial- vector fields 
(e.g., if a neutrino propagates in moving and polarized 
matter) can be treated with help of the classical field 
theory approach. 



V. CONCLUSION 

In conclusion we mention that the evolution of cou- 
pled classical fermions under the influence of external 
axial-vector fields has been studied in the present pa- 
per. We have discussed the particular case of two coupled 
fermions and formulated the Cauchy problem for this sys- 
tem. If the initial conditions were chosen in the appro- 
priate way, as it has been shown in Sec.[nJ the described 
system might serve as a theoretical model of neutrino fla- 
vor oscillations in matter. The initial conditions problem 
has been solved with help of the perturbation theory. We 
have found the zero and the first order terms in the fields 
distributions expansions over the external fields strength. 
It should be noted that the obtained results exactly take 
into account the Lorentz invariance as well as they are 
valid in (3+l)-dimensional space-time. The intensity of 
the zero order term corresponds to the case of the neu- 
trino flavor oscillations in vacuum. Therefore we have 
generalized our previous calculations performed in (1+1) 
dimensions in Ref. The first order correction is re- 
sponsible for the neutrino interaction with moving and 
polarized matter. Wc have obtained this intensity of the 
fermion field at great oscillations frequencies of the ini- 
tial field distribution, that corresponds to ultrarelativis- 
tic neutrinos. Note that we have compared our results 
with the transition probability formula for neutrino fla- 
vor oscillation in moving and polarized matter derived in 

Refs. mm and revealed an agreement in the case of 

weak external fields. This comparison proves the validity 
of the method elaborated in our work. Finally it has been 
demonstrated that neutrino flavor oscillations in moving 
and polarized matter could be described with help of the 
classical field theory. 

It is interesting to notice that along with the usual neu- 
trino oscillations phase equal to Am 2 /(4o;) the classical 
field theory approach also reveals rapid harmonic oscil- 
lations on the frequency [see, e.g., Eqs. (|3.4I) and 14.41) ] 



^rapid 



fiH + ^H 



4u 



However these terms are suppressed by the ratios mk/to 
which are small for great values of u). This case 
corresponds to ultrarelativistic neutrinos. The analo- 
gous terms were discussed in many previous publica- 
tions devoted to neutrino flavor oscillations (see, e.g., 
Refs. [6, 8]). In our work it has been demonstrated that 
such contributions to the neutrino transition probability 
appear even in the classical field theory approach. These 
terms arise from the accurate account of the Lorentz in- 
variance. 
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